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Observation of the Brownian motion of a small probe interacting with its en- 
vironment is one of the main strategies to characterize soft matter^"''. Essen- 
tially two counteracting forces govern the motion of the Brownian particle. 
First, the particle is driven by the rapid collisions with the surrounding sol- 
vent molecules, referred to as thermal noise. Second, the friction between 
the particle and the viscous solvent damps its motion. Conventionally, the 
thermal force is assumed to be random and characterized by a white noise 
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spectrum. Friction is assumed to be given by the Stokes drag, implying that 
motion is overdamped. However, as the particle receives momentum from the 
fluctuating fluid molecules, it also displaces the fluid in its immediate vicin- 
ity. The entrained fluid acts back on the sphere and gives rise to long-range 
correlation^'^. This hydrodynamic memory translates to thermal forces, which 
display a coloured noise spectrum. Even 100 years after Perrins pioneering 
experiments on Brownian motion^, direct experimental observation of this 
colour has remained elusive*. Here, we measure the spectrum of thermal 
noise by confining the Brownian fluctuations of a microsphere by a strong op- 
tical trap. We show that due to hydrodynamic correlations the power spec- 
tral density of the spheres positional fluctuations exhibits a resonant peak in 
strong contrast to overdamped systems. Furthermore, we demonstrate that 
peak ampliflcation can be achieved through parametric excitation. In analogy 
to Microcantilever-based sensors'''^" our results demonstrate that the particle- 
fluid-trap system can be considered as a nanomechanical resonator, where the 
intrinsic hydrodynamic backflow enhances resonance. Therefore, instead of 
being a disturbance, details in thermal noise can be exploited for the devel- 
opment of new types of sensors and particle-based assays for lab-on-a-chip 
applications^^'^^. 

Einsteins theory of Brownian motion published in 1905^^ received considerable attention, 
and was later reformulated in terms of a Langevin equation^"*. There, particle motion is driven 
by thermal fluctuations induced through collisions with the fluid molecules. These rapid kicks 
are assumed to be random and independent at frequencies much smaller than the collision rate 
of ~ 1 THz. The thermal force spectrum is then constant over a wide range of frequencies, 
which is commonly referred to as white noise spectrum^"^. However, when the densities of 
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particle and fluid are comparable, their coupling becomes important^' °. As the suspended 
particle fluctuates through the solvent, long-range correlations build up due to momentum ex- 
change, leading to hydrodynamic memory in the solvent by the particle's past motion. By the 
fluctuation-dissipation theorem, the statistics of the thermal force is characterized through a 
delta-correlated white-noise term, and a coloured frequency-dependent component, which re- 
flects the retarded viscous response of the fluid continuum to the particle. 

To directly measure the predicted correlations in thermal noise, we combined strong optical 
trapping with high-resolution interferometric detection^^. Under such conditions, the force bal- 
ance for the sphere, ignoring its inertia, reads -Ffr(t) — Kx{t) + -Fth(t) = 0, where x{t) is the 
bead's displacement from the trap center, K the stiffness of the harmonic optical trap, and Ff^{t) 
the friction force on the particle. At low frequencies, respectively long times, the trapping force 
dominates over friction, and the Langevin equation reduces to Kx{t) ~ -Fth(^)- Consequently, 
by monitoring the fluctuating motion of the bead, we effectively probe the thermal force of the 
fluid ^. At long times, correlations in thermal noise 

{F,i,{t)F,i,{0))c^K^x{t)x{0)), (1) 

become directly accessible by calculating the positional autocorrelation function, PAF(i) — 
{x{t)x{0)) from the recorded position fluctuations of the trapped sphere (Fig. la). Similarly, 
deviations from white noise in the power spectral density (PSD) of the displacement reflect the 
colour of thermal force ^. 

An overview of the characteristic time scales for trapped Brownian motion is displayed in 
Fig. lb. Comparing Stokes friction with the harmonic restoring force suggests the trap relax- 
ation time Tk = j/K as the longest characteristic time scale, where 7 = 6-KrjR is the stationary 
friction coefficient of the bead for macroscopic no-slip boundary conditions, rj the shear viscos- 
ity of the fluid and it! the radius of the sphere. Momentum is transferred from the particle to the 



3 



fluid at times Tp — m/j, with m the mass of the bead. Based on these contributions, and for not 
too strong trapping, tk > 4rp, the bead behaves like a harmonic oscillator in the overdamped 
regime. The corresponding PSD matches a Lorentzian, PSD (/) = 2(/i;BT/i^)rK/[l+(27r/rK)^] 
for frequencies around the comer frequency, /k = l/27rrK. Furthermore, the PAF displays an 
exponential relaxation PAF(t) = (fc^T/i^) cxp(— ^/tk) at time scales t ~ Tr. However, in a 
fluid this picture is incomplete, as hydrodynamic backflow has to build up first by vorticity dif- 
fusion (Fig. lb, middle), and the additional time scale, Tf — B?/ u, becomes important, where u 
denotes the kinematic viscosity of the fluid. The sequential order and magnitude of the various 
time scales depends on the size and mass of the Brownian particle, the nature of the solvent, 
and the stiffness of the optical trap. 

For a given solvent and particle, it is possible to investigate the dynamics of the system at 
different times by increasing K, thus shifting tk towards Tf or Tp. The experimental challenge 
consists in exploring a wide range of time scales. Therefore, we optimized our set-up^^ to 
obtain high trapping strength and best signal-to-noise ratio with a bandwidth covering at least 
four decades. The achieved resolution was on the order of Inm in space and close to l^s in 
time. 

In typical optical tweezers experiments, silica or polystyrene spheres immersed in water and 
with sizes < l|im are used. This yields values for Tf and Tp below l^s, thus, below our temporal 
resolution limit^^. Instead, we used melamine resin beads with a density pp — 1510 kg/m^ and 
a radius of 1 .0 ± 0. l|.im or 1 .5 ± 0. l|^m suspended in an acetone solution with u = 5 - 10~^m^/ s 
at 22°C. Then, Tp 1 i^s and Tf !^ 3 |is. Furthermore, the difference between the refractive 
indices of resin (n — 1.68) and acetone (n — 1.36) was high enough to provide a good trapping 
efficiency^*^ and allowed for the visualization of the particle by optical microscopy, which was 
not the case for e.g. silica particles (n = 1.37). With such an experimental setup, the magnitude 
of the trap relaxation time tk could be brought close to Tf, but still with tk > 4Tp. 



4 



The position fluctuations x{t) of a single resin sphere confined by the stiffest trap achievable 
with our set-up were acquired during T ~ 50 s at a sampling rate of 1 MHz. A comparison of 
PSD(/) at different experimental conditions (Fig. 2) confirmed that the choice of the bead size, 
fluid, and trap stiffness, hence Tf, respectively tk, determined the emergence of a resonance. 
When using water (red symbols), the maximal comer frequency /k we obtained resulted in 
an enhancement of the PSD close to our noise limit. Yet, striking deviations from a simple 
Lorentzian (red line), typical of overdamped systems, are clearly visible. The measured decay 
is steeper and shifted to higher frequencies. Replacing water by acetone decreased the ratio 
Tk/ Tf . As a consequence, a resonance peak appears indicating an enhanced spectrum of the 
thermal force with increasing frequencies. Thus, deviations from white noise are towards the 
blue at frequencies, which are by orders of magnitude smaller than the collision rate of the 
solvent molecules. In the time domain, the corresponding PAF(t) exhibits a zero crossing (Fig. 
3, green and blue symbols) resulting from anticorrelations, which are in remarkable contrast 
to the purely relaxing behavior of overdamped motion. This negative overshoot reflects the 
resonant peak in the PSD. The long-time decay is significantly slower than exponential and 
corroborates persistent correlations. Deviation from an exponential relaxation (green and blue 
lines) are substantial already at short times, characterized by a slower and then more rapid 
decay. 

For a quantitative description, we use a Langevin equation with no-slip boundary conditions 
accounting for slow vortex diffusion^^'^^ and trapping^ Our data for the PSD are in excellent 
agreement with the theoretical expression (black lines in Fig. 2), and were fitted simultaneously 
in the two dimensions lateral to the optical axis using the same Tf. We achieved trap stiffnesses 
of up to 345 |iN/m. The shape of the PSD is only determined by the ratio tk/ Tf , and hence 
reflects, at a given bead radius, the fluid properties through Tf , but is independent of the particle's 
mass. The fitted parameters were then used to compare the measured PAF(f) with theory^^. The 
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theoretical curves (black lines in Fig. 3) capture the data quantitatively over two decades in time 
and three orders of magnitude in signal, highlighting our unprecedented signal to noise ratio. 
We observe a hydrodynamic power-law tail PAF(i) ~ —{kBT^/K'^)-\jTi/A'Kt^, which by Eq. 
(1) directly reflects the corresponding persistent correlations in the thermal forces 

(i^th(t)Fth(O)) = -kBT-^^/^^ t-^/^ for i > . (2) 

As just shown, experimental access to short times reveals a resonance in Brownian mo- 
tion, where overdamped motion is commonly assumed. The main technical achievement is to 
strongly confine the particle and to detect its position fluctuations at highest bandwidth*^. For 
applications, it is desirable to enhance this effect. Stronger and narrower resonances can be ob- 
tained in the inertial regime, where Brownian motion is also sensitive to the particle's mass^^'^. 
To reach this window, Tk has to be approached to Tp by increasing K or m. Although heavier 
particles, which simultaneously allow for more efficient trapping are still to be developed^^, the 
timeline displayed in Fig. lb can be explored theoretically and by means of computer simula- 
tions. The appropriate Langevin equation, mx{t) = Ffr{t) — Kx{t) + Fth(t), comprises then 
a mass term, which eventually becomes larger than friction. By decreasing tk, the resonance 
is increased (Fig. 4a, continuous lines). Interestingly, compared to the simple overdamped 
harmonic oscillator (dashed lines), the peak is significantly enhanced and its position, /max. is 
shifted to lower frequencies by the contribution of hydrodynamic memory. 

The study of Brownian motion at short times in a medium displaying hydrodynamic effects 
has also become accessible by advanced simulation techniques. We employed multi-particle 
collision dynamics^^ with a molecular dynamics coupling between the solute and solvent par- 
ticles. The method yielded a highly compressible solvent and reflected correctly the hydrody- 
namic effects at coarse-grained scales. The approach was implemented most conveniently for 
full-slip boundary conditions at the solute-solvent interface, whereas our experiments obeyed 
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no-slip boundary conditions. The results of the PSD corresponding as closely as possible to our 
macroscopic system show that a resonance emerged in the friction-dominated, as well as in the 
inertial regime (Fig. 4b, full circles). However, weak coupling between the Brownian particle 
and the surrounding fluid yielded a much weaker resonance as opposed to no-sUp conditions 
(Fig. 4a). The collected data are in semi-quantitative agreement with the theoretical curves 
(Fig. 4b, coloured lines), where the friction coefficient is evaluated for a compressible fluid 
with full-slip boundary conditions'^. 

Further peak amplification can be achieved via parametric resonance'"*''^ by periodically 
modulating the trap strength, K{t) — K[l + cos(27r/exc^)]- We obtained the theoretical 
PSDexc(/) normalized to the initial value of the non-excited PSD from the solution of a para- 
metrically modulated Langevin equation, including hydrodynamic memory, by second order 
perturbation theory in the reduced modulation amplitude g. Similar to the harmonic oscillator, 
also in the presence of coloured friction best excitation is achieved at a frequency /exc — S/max, 
yielding an increase of up to 20% when g = 50% (Fig. 4c, empty circles). Comparable results 
were obtained within computer simulations (empty circles in Fig. 4b). 

In view to develop particle-based assays^ ^'^', we evaluated the sensitivity of the peak by 
slightly changing the size of the particle. As shown in the inset of Fig. 4c, a change of 5% in 
the bead radius results in a resonance frequency shift of 1 kHz. We anticipate that changes in the 
particle's morphology, like swelling or a reaction occurring at the surface of such a Brownian 
probe will alter its short-time dynamics and become detectable. As also single cells, microor- 
ganisms and microcarriers can be bound harmonically short-time detection of their Brow- 
nian fluctuations may become a sensitive way to characterize their state or evolution in native 
solutions and without specific markers. Reciprocally, also changes in the medium surrounding 
the probing particle modulate its fluctuation spectrum'"'*, and a dynamic polymer system may 
be studied in great detail. 
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METHODS SUMMARY 

Melamine resin spheres (Sigma-Aldrich, Analytix 5, 2001) were suspended in high purity 
acetone at minimal concentrations to allow trapping and observation of exclusively one single 
particle. The sample chamber was made of a microscope slide having a ground-in polished 
cavity in its center (18mm in diameter and 0.8mm in depth) (Marienfeld GmbH, Germany). 
The cavity was filled with the bead-acetone mixture and sealed with a glass coverslip (thickness 
~ 130|am) and high viscosity vacuum grease to prevent evaporation. Special care was taken not 
to introduce bubbles into the sample chamber. After loading, the sample chamber was mounted 
onto our custom-made inverted microscope-optical trap set-up. 

Since the trap was asymmetric^^, we obtained different Tk for each dimension from fitting 
our data to the theory. Both values differed by maximally 20%. The determined values for T{ and 
Tk were used to calculate u and K. The bead radii were chosen to be within the tolerances given 
by the manufacturer. The slightly higher acetone viscosities of i/ = 5 • 10~^m^/s at the labora- 
tory temperature of 22°C, compared to the literature value ofu = A - 10~^m^/s were attributed 
to a contamination of our acetone stock solution by atmospheric water during experimentation. 
In situ viscometry measurements^^ were performed to cross-check for the value of the acetone 
viscosity obtained through the determination of Tf. We obtained v = (4.8 ± 0.6) • 10~'^m^/s for 
10 measurements with different beads. As acetone has a very low absorption coefficient at our 
laser wavelength A = 1064 nm, no heating of the solution was expected. 

Full Methods and any associated references are available at the end of the manuscript. 
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FIGURES AND LEGENDS 
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Figure 1 | Characteristic time scales for a Brownian particle confined by the harmonic 
potential of an optical trap, a, Trajectories of the trapped sphere measured in the 2 dimen- 
sions lateral to the optical axis, at three different time intervals, b, At very short times t <^ 
the particle undergoes ballistic motion (left). At the hydrodynamic time scale Tf hydrodynamic 
backflow develops (center; solid lines show the emerging fluid velocity field, arrows are ob- 
tained as part of the present computer simulations). Finally, for t > the harmonic potential 
of the optical trap sets in and confines the diffusion of the particle (right). 
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Figure 2 | Emergence of a coloured peak in the thermal noise spectrum. Log-linear plot 
of the PSD normalized to its zero-frequency value 2kBTr}^/ K of a melamine resin sphere in 
pure water (red: T{ = 2.1|as, tk = 130|as, R = 1.42 |am, u = 9 ■ 10"W/s, k = 196|jN/m) 
or in acetone solution (z/ = 5 ■ 10~^m^/s) with increasing Tf/r^ ratios (magenta: Tf = 2.3|as, 
Tk = 75.9|as, R = 0.96 |am, K = 77|jN/m; green: rf = 5|as, tk = 45|as, R = 1.42|am, 
K = 190 |aN/m; blue: Tf = 2.4|as, tk = 19.2|as, R = 0.99 |am, K = 309 |aN/m). Data were col- 
lected for T ~ 50 s, corresponding to ~ 5 • 10^ datapoints. Error bars indicate the standard error 
on the mean. The PSD were computed from overlapping windows of 2^^ points, blocked in 100 
bins per decade. The black lines correspond to the hydrodynamic theory^^ The parameters Tf 
and Tk were extracted from the fit. The coloured lines are Lorentzian spectra corresponding to 
non-hydrodynamic systems. Inset: Zoom into the enhancement of the PSD, which reflects the 
colour of thermal noise (rK/rf = 8, blocked in 20 bins per decade). 



14 




Time t[s] 

Figure 3 | Anticorrelations in tlie PAF reveal tliat Brownian motion is not overdamped. 

Double logarithmic representation of the normalized PAF(t) for two datasets of Fig. 2 (blue: 
TK/Tf = 8, green: r^/rf = 9). The black lines correspond to the full hydrodynamic theory 
including inertial effects^^. The persistent anticorrelations are visible after the zero crossing 
(narrow spike) and approach a t^^/^ power law decay. The blue, respectively green lines indi- 
cate exponential relaxations and serve as guide to the eye. Inset: Semilogarithmic plot of the 
same data. 
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Figure 4 | Amplification of the resonant peak, a, Theoretical PSD^' for a resin sphere in ace- 
tone {R = 1.5|am, Tp = 2.3|js, Tf = 5.4|as) for very strong traps (K = 1.02 , 2.04 , 4.07 mN/m; 
blue, green, and red graphs, respectively). The dashed lines correspond to the damped harmonic 
oscillator Inset; Corresponding plot of the fluid's velocity field developing around the sphere 
moving along the x-direction for no-slip boundary conditions at / = l/(27rrf). The arrows in- 
dicate the direction of the velocity field, while the colours denote the velocity magnitude rising 
from green to red. b, Simulation data (full circles) for an equivalent system with Tp = 0.67rf 
evaluated in a compressible fluid for full-slip boundary conditions, 7 = Aivqa. The coloured 
theoretical lines account for vortex diffusion as well as for sound waves^^. The empty circles 
are simulation data for an excitation with g = 50% and /exc = 2/niax- Inset; Corresponding 
velocity field at / = l/(27rrf). c, PSDexc(/) for the same particle as in a with g = 50% and 
/exc — 2/max (empty circles). Inset; Shift of the peak frequency when increasing the radius of 
the particle by 5% respectively 75nm. All curves are normalized to the zero-frequency value of 
the respective non-excited PSD. 
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FULL METHODS 

Instrument Set-up. Our Instrument consists of a custom-built inverted light microscope 
with a three-dimensional (3D) sample positioning stage, an infrared laser for optical trapping, 
and a quadrant photodiode for high-resolution and time-resolved position detection. The trap- 
ping beam is produced by a diode-pumped, ultra-low noise Nd:YAG laser with a wavelength 
of A = 1064nm (IRCL-500-1064-S, CrystaLaser, USA), and a maximal light power of 500mW 
in continuous wave mode. Best trapping efficiency is achieved by expanding the effective laser 
beam diameter 10 times with a telecentric lens system (Sill Optics, Germany). Further details 
on the set-up are given in^^ and its online supplement. 

Data analysis. Since the location of the trap center is not known with high precision, the 
position autocorrelation function, PAF(t) = {x(t)x{0)) cannot be computed directly. Rather 
the RAF is derived from the mean-square displacement, MSD(t), via the relation PAF(t) = 
ksT/K - MSD(i) /2. The MSD is computed from the position time trace x{t) as 



with j = 0, N,tj = jAt, and fs = l/At= IMHz our sampling frequency. Using the long- 
time limit of the MSD determined by the confinement of the trap, MSD(t ^ oo) = 2kBT/K, 
the normalized RAF reads 



We decomposed x{t) into Fourier modes yielding the Fourier transform 

rr/2 

xtU)= / x{t)ex^{2nift)dt, (5) 

J-T/2 

where the frequencies / are integer multiples of 1/T. The power spectral density is then ob- 
tained for long observation times as 




(3) 



R\F(i)/R\F(0) = 1 - MSD(i)/MSD(oo). 



(4) 




(6) 
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When plotting the PAF and the PSD on a logarithmic time, respectively, frequency scale, further 
data processing can be performed to extract smooth functions from noisy data^\ The data are 
collected in equidistant time steps and since both PSD(/) and PAF(t) vary essentially only on 
a logarithmic horizontal axis, the number of points in a log-plot will increase for long times or 
low frequency, respectively. Therefore, data are commonly averaged by 'blocking' consecutive 
data points, resulting in equidistantly distributed points on the logarithmic scale. The data 
scattering around their average value gives the standard error on the mean. Blocking allows to 
fit the data by least-square methods, which assumes that the analyzed data points are statistically 
independent and conform to a Gaussian distribution. 

Simulations. The simulations are based on the method of multi-particle collision dynamics 
(MPC)^^'^^. This mesoscale simulation technique models the solvent as a set of point particles 
of mass mf and density pf with continuous coordinates r and velocities v. At every simulation 
step the solvent particles undergo propagation and collision. In the propagation step, the coor- 
dinates of each particle are updated according to f{t + 5t) = f{t) + v{t)5t. The particles are 
then sorted into collision cells, and their velocities are rotated with respect to the center-of-mass 
velocity of the respective collision cell by a fixed angle -& around a random axis. The direction 
of the rotation axis is the same within each collision cell, but is statistically independent be- 
tween different cells and for any cell at different times. The standard grid-shifting procedure 
to guarantee Galilean invariance is employed. Mass is measured in units of the solvent particle 
mass ruf, length in units of the collision cell size b, energy in units of the thermal energy /c^T 
and time in units of to = ^/nifb'^/kBT. Our model consists of a single colloidal particle of mass 
m — 335mf embedded in a sea of solvent particles with density pf — lOnif/b^ within a simu- 
lation box with a lateral size L = 32b with periodic boundary conditions. The solute-solvent 
particle coupling is given by a truncated Lennard- Jones potentiaP^ with the interaction radius 
a — 2.0b. Since no angular momentum is transferred to the particle by the collisions, we simu- 
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late effectively full-slip boundary conditions. In order to reproduce fluid-like behavior, we use 
the collision time step O.OSto and collision angle d = 130° The kinematic viscosity is then 
calculated^^ XobQ v — 77/pf = and the Schmidt number Sc — v/D = 65. The dif- 

fusion coefficient D of the colloidal particle measured simultaneously through its mean-square 
displacement and the integral of the velocity auto-correlation function is D = 0. 003 (2)6^ /to- 
Then, the effective hydrodynamic radius was determined io R = 1.596 by the Stokes-Einstein 
relation D — ksTl^, where 7 = 47777 for slip. Then the vortex diffusion time scale is 
calculated to Tf = 1.51io, whereas the momentum relaxation time scale Tp = m/7 evaluates 
to Tp = l.Oto- The colloidal diffusion time td = d} jD = MHq is much larger in agreement 
with the experimental conditions. Due to the high compressibility of the solvent, the sonic time 
Ts = R/cq = 1.59to, with the isothermal sound velocity cq = ^JksT jmi, becomes comparable 
to vortex diffusion. Ignoring the sonic scale, the resulting order of time scales follows the right 
physical hierarchy Tp < Tf <^ td present in the experimental system, as prescribed in^^. The 
force exerted on the colloidal particle by the optical trap is taken to be of a harmonic form, 
-ftrap(^) — —Kx{t), where K is the spring constant of the trap and x{t) is the displacement of 
the particle relative to the center of the trap. The trapping potential gives rise to an additional 
characteristic time scale = ^/K. In the simulations for parametric excitation, the spring 
constant of the trap is modulated as K{t) = K[l + g cos(27r/exc^)], with a relative modulation 
amplitude g — 0.5. To prevent heating, the system is thermalized with a time step itherm = 
by rescaling the velocities of all particles such that the temperature of the system corresponds 
to the thermal energy ksT. 

Theory. The Langevin equation for a particle subject to Brownian motion in a harmonic 
potentiaF^ reads 

mx{t) = Ffr(t) - Kx{t) + Fth(t), (7) 
where x{t) is the displacement of the bead relative to the trap center. The mass of the particle 
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is again denoted by m, K represents the stiffness of the trap, Ffr{t) the deterministic friction 
forces, and Fih{t) is the thermal noise. Since the hydrodynamic flow pattern has to build up 
first, the friction force is retarded 

Ffrit)^- J -f{t-t')x{t')dt', (8) 

where the kernel ^{t) vanishes for i < by causality. By the well-known fluctuation dissipation 
theorem^^, the time correlation function of the random forces is determined by the same friction 
kernel 

{F^{t)F^{t'))^kBT^{\t-t'\), (9) 

and ksT is the thermal energy. It is convenient to perform a temporal Fourier transform, con- 
vention x{uj) — J e^'^*x{t)dt, where u — 27r/ denotes the angular frequency. The Langevin 
equation is then equivalent to 

[-muj^ - iuj^{uj) + K]x{uj) = F^{uj). (10) 

For the correlation function of the random forces in the frequency domain one obtains 

(Fth(c^)Fth(a;0*) = A7rkBTRc[^uj)]5iu - ou'). (11) 

For finite observation times T, one replaces 27rd{u} — u') T, and infers the power spectral 
density of the forces as 2kBTRe[y{u!)]. The colour of Brownian motion becomes manifest in 
the frequency dependence of the friction kernel. Defining the Green function G{uj) = [—mw'^ — 
+ the power spectral density (PSD) corresponding to the positional fluctuations 

(PAF) of the bead is then calculated to 

PSD(a;) = 2A;BTRe[7(a;)]|G(a;)|2 = 2kBT^^^^^. (12) 

(jO 

By the Wiener-Khinchin theorem the PAF(t) = {x{t)x{Q)) is then obtained from the corre- 
sponding PSD by an inverse Fourier transform 

PAF(t) = /— e-^'^*PSD(a;). (13) 
J 27r 
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Einstein's theory is recovered by assuming a white noise behavior at the frequencies of in- 
terest for the power spectral density of the forces 2kBTRe[^{uj)] ^ 2/0^X7, with 7 = 7(0; = 0). 
This corresponds to a time correlation function that is short lived (Ffr(i)Ffr(f)) — 2kBT^5{t — 
t') as is used in many introductory textbooks. For a spherical particle and no-slip boundary 
conditions the Stokes friction coefficient equals 7 = 6nr]R, where R is the radius of the bead. 
The assumption that the friction force is given by Stokes' law is incorrect, since the long-range 
flow profile has to build up first. The frequency-dependent friction coefficient 7(0;) correspond- 
ing to a sphere performing small-amplitude oscillations in an incompressible fluid with no-slip 
boundary conditions was already calculated by Stokes 



Here a — and the square root is to be taken with positive real part, whereas Tf = R? jv 

denotes the vortex diffusion time. Since 67rr]RTf/9 — mf/2, the last term can be interpreted as 
effectively increasing the mass of the particle by half the mass, mf, of the displaced fluid. The 
square root describes the evolution of the backflow by vortex diffusion and is non-analytic in 
frequency. In multi-particle collision dynamics, the large compressibility of the simulated fluid 
gives rise to an additional time scale = R/cq, with cq the isothermal sound velocity, which 
comes close to the vortex diffusion time Tf. Furthermore, the simulations employ effectively 
full-slip boundary conditions, since no angular momentum is transferred to bead during the col- 
lisions with the solvent particles. The corresponding frequency-dependent friction coefficient 
was calculated in^'^'^^^? 



Here X = R/ -^Cq — iu!{( + 4rj/3)/pf accounts for the sound propagation, with ( the bulk vis- 
cosity and pf the mass density of the fluid. One readily checks that 7(0; = 0) = AnrjR assumes 



7(0;) = 67ir]a{l + a + 



(14) 



ARnr] (1 + A) (18 + 18a + Sa^ + + 4(1 + a)X'^ 
~3 (2 + 2A + A2)(3 + a) + 2(1 + a)X'^/a'^ 



(15) 
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the correct value for the stationary friction of a sphere for full slip. 
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